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Abstract. In the Hartree approximation of Cornwall- Jackiw-Tomboulis (CJT) 
formalism of the real scalar field theory, we show that for the strongly coupled scalar 
system near phase transition, the shear viscosity over entropy density is small, however, 
the bulk viscosity over entropy density is large. The large bulk viscosity is related to 
the highly nonconformal equation of state. It is found that the square of the sound 
velocity near phase transition is much smaller than the conformal value 1/3, and the 
trace anomaly at phase transition deviates far away from 0. These results agree well 
with the lattice results of the complex QCD system near phase transition. 



Strongly coupled matter near phase transition 



2 



It has been beheved that quark matter created at RHIC behaves hke a nearly 
"perfect" fluid. One crucial quantity to identify this property is the small ratio of 
shear viscosity over entropy density t]/s. In order to fit the elliptic flow at RHIC, the 
hydrodynamic simulation shows that a very small shear viscosity is required [Ij . Lattice 
QCD calculation confirmed that rj/s for the purely gluonic plasma is rather small and 
in the range of 0.1 — 0.2 |2J. 

However, recent lattice QCD results showed that another important transport 
coefficient, the bulk viscosity over entropy density ratio (/s rises dramatically up to 
the order of 1.0 near the critical temperature [Hll]. The sharp peak of bulk viscosity 
at Tc has also been observed in the linear sigma model [5], and the increasing tendency 
of C/s below Tc has been shown in a massless pion gas [6]. The large bulk viscosity near 
phase transition is related to the nonconformal equation of state [7J. The correlation 
between large bulk viscosity and the nonconformal euqation of state has been shown in 
Ref. 0. 

The bulk viscosity is related to the correlation function of the trace of the energy- 
momentum tensor 6f,: 

C = -lim-/ dt d're^-'{[9^X^),9^^m). (1) 

y t<;-*0 to Jo J t I- 

According to the result derived from low energy theorem, in the low frequency region, 
the bulk viscosity takes the form of [3] 

C = 7T^ \t'^ ^"^ ~J/^^ + 16|6.|| = {-16e + 9Ts + TC.} . (2) 



dT ' "'J 9cuo 

Here e is the total energy density of the system, e„ the negative vacuum energy density, 
s is the entropy density and is the specific heat. We have introduced the normalized 
pressure density pt and energy density e^. The parameter ujq = uJo{T) is a scale at 
which the perturbation theory becomes valid. 

There are several quantities to characterize the confomality of the system. One is 
the square of the speed of sound cj., which is related to Pt/^t and has the form of 

2 ^ dp ^ S ^ ^ 

' de Tds/dT C ^ ^ 

At the critical temperature, the entropy density as well as energy density change most 
fastly with temperature, thus one expect that cj. should have a minimum at T^. Another 
quantity to describe the conformality of the system is the trace anomaly of the energy- 
momentum tensor Tf^^, i.e, 

^ = J^^-^J^ = T^iPT/T% (4) 

We define A/rf as the "interaction measure", with d the degeneracy factor. 

The conformal limit has attracted much attention in recent years, since people 
are trying to understand strongly interacting quark-gluon plasma by using AdS/CFT 
techniques. In conformal field theories including free field theory, pt/^t = = 1/3, 
A = 0, and the bulk viscosity ( is always zero. Lattice results show that at 
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asymptotically high temperature, the hot quark-gluon system is close to a conformal 
and free ideal gas. 

However, lattice results show that near deconfinement phase transition, the hot 
quark-gluon system deviates far away from conformality Both Pt/^t and 

show a minimum around 0.07, which is much smaller than 1/3. For the SU{3) 
pure gluon system, the peak value of the trace anomaly A^^-^ reads 3 ~ 4 at T^ax 
and the corresponding "interaction measure" is 

A^Ar/dc = 0.2 ~ 0.25, with the 
gluon degeneracy factor dc = 16. (Note that here Tmax — l-lT'c is the temperature 
corresponding to the sharp peak of A.) For the two- flavor case, the lattice result of 
the peak value of the trace anomaly A^^^^ reads 8 ~ 11, the corresponding interaction 
measure at T^ax is given as A^^t^ /{do + dq) = 0.28 ~ 0.4, with quark degeneracy 
factor dq = 12. 

Due to the complexity of QCD in the regime of strong coupling, results on hot quark 
matter from lattice calculation and hydrodynamic simulation are still lack of analytic 
understanding. In recent years, the anti-de Sitter /conformal field theory (AdS/CFT) 
correspondence has generated enormous interest in using thermal A/" = 4 super- Yang- 
Mills theory (SYM) to understand sQGP. The shear viscosity to entropy density ratio 
rj/ s is as small as l/47r in the strongly coupled SYM plasma [9]. However, a conspicuous 
shortcoming of this approach is the conformality of SYM: the square of the speed of 
sound cj. always equals to 1/3 and the bulk viscosity is always zero at all temperatures 
in this theory. Though (^/s at is non-zero for a class of black hole solutions resembling 
the equation of state of QCD, the magnitude is less than 0.1 [10], which is too small 
comparing with lattice QCD results. 

It has been found in Ref. [H] that in the simplest real scalar model with Z{2) 
symmetry breaking in the vacuum, rj/s behaves the same way as that in systems of 
water, helium and nitrogen in first-, second-order phase transitions and crossover |12j . 
In Ref. |T3], we have investigated the equation of state and bulk viscosity in the real 
scalar model, and compare the result in this simplest relativistic system with that of 
the complex QCD system. 

The Lagrangian of the real scalar field theory has the form of 

C = l{d,<l>r - \act>' - \b<p\ (5) 

with a the mass square term and b the interaction strength. This theory is invariant 
under —(f) and has a Z2 symmetry. In the case of a < and b > 0, the vacuum at 
T = breaks the Z2 symmetry spontaneously. The Z{2) symmetry will be restored at 
finite temperature with a second-order phase transition. 

At finite temperature, the naive perturbative expansion in powers of the coupling 
constant breaks down. A convenient resummation method is provided by the extension 
of Cornwall- Jackiw-Tomboulis (CJT) formalism [Tl] to finite temperature. The CJT 
formalism is equivalent to the ^-functional approach of Luttinger and Ward ^5J and 
Baym [16]. In our calculation, we only perform the Hartree approximation for the 
effective potential, i.e, only resum tadpole diagrams self-consistently and neglect the 
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exchange diagrams. The effective potential in the CJT formalism reads [T7j 

mS] = ^l^[\nS-\K) + S^\K)SiK)-l]+ V^^, S] + Ui^) , (6) 

where U{(f)) = a/2 0^ + 6/4 0^ is the tree-level potential, and the 2PI potential 
V2[0, 5] = S{K,(j))^ in the Hartree approximation. S{So) is the full (tree-level) 

propagator and takes the form of S~^{K, 0) = — ii'^ + m^(0), Sq^{K, 0) = — K'^ + ml{(j)) 
with the tree-level mass = a + 36 0^. 

The gap equations for the condensation 0o and scalar mass m are determined by 
the self-consistent one- and two-point Green's functions 



6Q 



= ^ 

— ' 5S 

4>=<t>o,S=S{4>o) 



^ 0. (7) 

!.o,5=5(</.o) 

The entropy density is determined by taking the derivative of effective potential 
with respect to temperature, i.e, s = —dQ{(j)o)/dT. In the symmetry breaking 
case, the vacuum effective potential or the vacuum energy density is negative, i.e, 
Qy = f2(0o)|T=o < 0. The normalized energy density and pressure density pt can be 
calculated by px = —^t with Qt = ^(0o) ~ and = —pT + Ts. 

As a reference for complex QCD system, we investigate the equation of state and 
transport properties for the real scalar field theory with Z{2) symmetry breaking in the 
vacuum and 2nd order phase transition at finite temperature. The trace anomaly A the 
specific heat Cy as well as bulk viscosity to entropy density ratio (/s show upward cusp 
at Tc, and their peak values increase with the increase of coupling strength. The ratio of 
pressure density over energy density Pt/^t and the square of the sound velocity show 
downward cusp at T^, which is similar to the behavior of rj/s found in Ref. [11], and 
the cusp values decrease with the increase of coupling strength. These cusp behaviors 
at phase transition resemble lattice QCD results. In Fig. [1], we only show the trace 
anomaly A and the ratio of bulk viscosity to entropy density ratio (/s as functions of 
T/Tc for different coupling strength b. 
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Figure 1. The interaction measure [ct — ■^Pt)/T^ (^^ft) '^rid the bulk viscosity over 
entropy density ratio C./s (right) as a function of temperature T/Tc for different 
coupling strength b. 
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In the weak coupling case when b = 0.3, the cusp values of Pt/^t and at 
are close to the conformal value 1/3, both the trace anomaly A and the bulk viscosity 
to entropy density ratio (/s at Tc are close to conformal value 0. However, the shear 
viscosity over entropy density ratio r]/s is around 2000, which is huge comparing with 
the AdS/CFT limit l/4n. Here we have used the method in Ref.|llj to derive rj/s. To 
our surprise, we find that when b = 30, the strongly coupled scalar system can reproduce 
all thermodynamic and transport properties of hot quark-gluon system near T^. Pt/^t 
at Tc is close to the lattice QCD result 0.07, A/d = 0.48 {d = 1 for scalar system) at Tc 
is close to the lattice result of the peak value A^^^^/ {do + dq) ~ 0.4 at T^ax- The bulk 
viscosity to entropy density ratio C/s at Tc is around 0.5 ~ 2.0, which agrees well with 
the lattice result in Ref. (Note, here C,/s = 0.5,2 correspond to t^o = lOT, 2.5T, 
respectively.) More surprisingly, the shear viscosity over entropy density ratio r]/s at 
Tc is 0.146, which also beautifully agrees with lattice result 0.1 ~ 0.2 in Ref. [2]. In 
Table [H we compare our results of equation of state and transport properties in scalar 
field theory at Tc, and corresponding results in lattice QCD calculations [31 HJ [7], the 
Polyakov-loop Nambu-Jona-Lasinio (PNJL) model [HI [19], and black hole duals |10j . 
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Table 1. Thermodynamic and transport properties at T/Tc = 1 in scalar theory at 
weak coupUng b = 0.3 and strong couphng b = 30, in lattice QCD [3l HI [7], PNJL 
model [m \T9\ , and black hole dules [10] . The degeneracy factor d = 1 for real scalar 
model. 

In summary, in the Hartree approximation of CJT formalism, we have investigated 
the equation of state and transport properties of the real scalar field model with Z{2) 
symmetry breaking in the vacuum and 2nd-order phase transition at finite temperature. 
We have seen that at phase transition, the system either in weak coupling or strong 
coupling shows some common properties: 1) Pt/^t, the square of the speed of sound 
Cg as well a.s rj/s exhibit downward cusp behavior at Tc. 2) The trace anomaly A, the 
specific heat C„ as well as (/s show upward cusp behavior at Tc. The cusp behavior 
is related to the biggest change rate of entropy density at Tc. At weak coupling, the 
scalar system near phase transition is asymptotically conformal. However, the shear 
viscosity is huge. At strong coupling, the scalar system near phase transition is highly 
non- conformal, the shear viscosity is small, but the bulk viscosity is large. We can expect 
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that for any nonconfomal field tfieory, the corresponding system at strong coupling near 
phase transition exhibits highly nonconformality. AdS / CFT method maybe cannot help 
us understand the strongly interacting quark gluon plasma. Unexpectedly, the simplest 
scalar field model can do the job. We have found that lattice QCD results on the 
equation of state and transport properties near phase transition can be amazingly very 
well described by the simplest real scalar model at strong coupling when 6 = 30. It is 
urgent to include the bulk viscosity correction and the nonconformal equation of state in 
hydrodynamics to investigate hadronization and freeze-out processes of QGP created at 
heavy ion collisions [20l[2T], it would be also interesting to investigate how bulk viscosity 
affects charm radial flow [22] at RHIC. 
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